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ABSTRACT 

We investigate the flux and event rate of supernova relic neutrinos (SRNs) and discuss their implications for 
the cosmic star formation rate. Since SRNs are diffuse neutrino background emitted from past core-collapse 
supernova explosions, they contain fruitful information on the supernova rate in the past and present universe, 
as well as on the supernova neutrino spectrum itself. As reference models, we adopt the supernova rate model 
based on recent observations and the supernova neutrino spectrum numerically calculated by several groups. 
In the detection energy range > 10 MeV, which will possibly be a background-free region in the near 
future, the SRN event rate is found to be 1-2 yr~' at a water Cerenkov detector with a fiducial volume of 
22.5 kton, depending on the adopted neutrino spectrum. We also simulate the expected signal with one set of 
the reference models by using the Monte Carlo method and then analyze these pseudodata with several free 
parameters, obtaining the distribution of the best-fit values for them. In particular, we use a parameterization 
such that /?sn(z) = ^sn(^ "'"^)"' where /?sn(z) is the comoving supernova rate density at redshift z and R^j^ and a 
are free parameters, assuming that the supernova neutrino spectrum and luminosity are well understood by way 
of a future Galactic supernova neutrino burst or the future development of numerical supernova simulations. 
The obtained Itr errors for these two parameters are found to be 6a/ (a) = 30% (7.8%) and (^TJ^f^/ {R^j^) = 
28% (7.7%) for a detector with an effective volume of 22.5 kton x 5 yr (440 kton x 5 yr), where one of the 
parameters is fixed. On the other hand, if we fix neither of the values for these two parameters, the expected 
errors become rather large, 6a/ {a) = 37% and iJ^sn/(^sn) ~ 55%, even with an effective volume of 440 
kton X 5 yr 

Subject headings: diffuse radiation — neutrinos — galaxies: evlution — supernovae: general 



1. INTRODUCTION 

In recent years we have made remarkable progress in our 
knowledge concerning how the cosmic star formation his- 
tory proceeded in the past and concerning the fraction of 
baryons locked up in stars and gas in the local universe. 
These points were inferred from observations of the light 
emitted by s t ars of various masses at various wavelengths. 
iMadau et alJ Jl996h investigated the galaxy luminosity den- 
sity of rest-frame ultraviolet (UV) radiation up to z 4, 
and they converted it into the cosmic star formation rate 
(SFR). The rest-frame UV light is considered to be a di- 
rect tracer of star formation because it is mainly radiated 
by short-li ved m assive stars. After the pioneering study 
by Madau et al.L a wealth of data have become available in 
the form of the cosmic SFR in a wide range of redshifts; 
these data were inferred from observations using far-infrared 
(FIR)/ s ub millimeter dust emission (Hughes et al. 1998; 
iFlores et al. il999i) and near-infrared (NIR) Ha line emis- 
sion (Galle go et 31111 995: Gronwall 1998; Tresse & Maddox 
[1^98; Glazebrooket al. 1999), as well as the rest-frame UV 
emiss io n from massive stars (Lillv et al. 1996; Cowie et al. 
T996t ICoiinollv et alJ ll997t~1Sa wicki. Lin. & Yee 1997; 



Trever et al.l 11998^ 'Madau. ^ozzetti. & Dickinson H 
Pascarelle. Lanzetta. & Fernandez-Soto 1998; Steid el et al.l 
19991) . 

In these traditional methods, however, there are a fair num- 
ber of ambiguities whe n the actual observables are converte d 
into the cosmic SFR JSomerville. Primack. & Fa"be3l200l . 
First, observable samples are generally flux-Umited, and thus 
the intrinsic luminosity of the faintest objects in the sample 

Electronic address: |Mdo@utap.phys.s.u-tokyo.ac.jpl 



changes with redshift. In order to understand the true red- 
shift dependence of the total luminosity density, this incom- 
pleteness is generally corrected by using a functional form 
(i.e., a Schechter function) of the luminosity function ob- 
tained from the observations themselves. Unfortunately, since 
the luminosity function is not well established observation- 
ally (especially for high-z regions), it is uncertain whether the 
Schechter-function fit is good enough or not. Second, the con- 
version from luminosity density to SFR generally relies on 
stellar population models, an assumed star formation history, 
and an initial mass function (IMF), which are not also well es- 
tablished yet. Finally, if the tracer of star formation is an op- 
tical or UV luminosity, then the effects of dust extinction are 
nonnegligible. Although this problem is less critical in other 
wave bands such as NIR Ha or FIR/submillimeter, the bulk 
of current data consists of rest-frame UV observations, espe- 
cially of high-redshift regions. After adopting some correc- 
tion law for dust extinction, the rest-frame UV data become 
rather consistent with Ha or submillimeter data points; still, 
in this case it is unknown whether the UV and submillimeter 
sources are identical, which is very important for measuring 
the cosmic SFR. 

Thus, our knowledge concerning the cosmic SFR is quite 
crude, and therefore another type of observation that is inde- 
pendent of the above methods would be very important. In 
this paper we consider supernova relic neutrinos (SRNs), i.e., 
a diffuse background of neutrinos that were emitted from past 
supernova explosions. Type lb, Ic, and II supernova explo- 
sions are considered to have traced the cosmic SFR, because 
they are directly connected with the death of massive stars 
with M > 8M0 and their lifetime is expected to be very short 
compared with the time scale of star formation. These events 
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' in the energy range > 19.3 MeV 
This severe constraint is only about 



are triggered by gravitational collapse, and 99% of the grav- 
itational binding energy is released as neutrinos; this basic 
scenario was roughly confirmed by the well-known observa- 
tion of the neutrino burst fr om SN 1987A by the Kamiokan de 
II and 1MB detectors (Hirata et al. 1987; Bionta et al. 1987). 
The advantages of using SRNs to probe the cosmic SFR are 
as follows: First, as mentioned above, if we can probe the su- 
pernova rate from SRN observation, it can be directly trans- 
formed to the SFR assuming the IMF, because supernovae 
are short-lived astrophysical events. The second advantage, 
which is more important, is that neutrinos are completely free 
of dust extinction. This point is the same with observations in 
the submillimeter wave band; however, neutrinos are emitted 
directly from stars, whereas submillimeter radiation comes 
from dust and is an indirect process. 

The SRN flux and the event rates at a currently working 
large-volume water Cerenkov detector, Super-Kamiokande 
(SK), have been investigated by many researchers using the- 
oretically/obser vationally mod el ed S FRs (Totani & Sato 
1^95; Totani. Sato. & Yoshii 1996; Malanev 1997; 
Hartmann & Wooslev 1997; Kaclinghat. Steigman. & Walker 
2000; Ando, Sato, & Totani 2003). More recently, flie SK 
collaboration obtained a 90% CL upper limit on the SRN 
flux, i.e., 1.2 cm~^ s 
jMaleketal.ll2()0l . 
factor of 3-6 larger than the typical theoretical models and 
is very useful for obtai ning several rough estimations of 
the cosmic SFR (Fukugita & Kawasaki"2003t IStrigari et alJ 
|2P03) and probing the properties of neutrinos as elementary 
particles (Ando & Sato 2003a; Ando 2003). 

However, we need a further ^ 40 years to reduce the cur- 
rent limit by a factor of 3 if we use the SK detector with 
current performance. This is because there is no energy win- 
dow for SRN detection where the SRN signal dominates other 
background events coming from various sources, such as so- 
lar, reactor, and atmospheric neutrinos, as well as cosmic- 
ray muons jAndo et al. 2003). Therefore, current observa- 
tions are seriously affected by the other backgrounds and 
take much time to reach the required sensitivity. In order 
to overcome this difficulty, a very interesting and promis- 
ing method was proposed to directly tag electron antineu- 
trinos (Df), and it is now in th e research and development 
phase JBeacom & Vagin3l2003 h. The basic idea is to dis- 
solve 0.2% gadolinium trichloride (GdCh) into the pure wa- 
ter of SK. With this mixture, 90% of the neutrons produced 
by the D^p — > e'^n reaction are captured on Gd and then de- 
cay with 8 MeV gamma cascades. When we detect these 
gamma cascades, as well as the preceding Cerenkov radiation 
from positrons, it indicates that these signals come from orig- 
inal not from other flavor neutrinos or muons. With this 
method, we can remove the background signals in the energy 
range 10-30 MeV, in which before removal there is a huge 
amount of background from solar neutrinos (I'g) and atmo- 
spheric muon-neutrinos (f^, P^j} or cosmic-ray muon induced 
events. Because the expected SRN rate is estimated to be 1-2 
yr"' in the energy range 10-30 MeV, the Gd-loaded SK detec- 
tor (Gd-SK) would enable us to detect a few SRN events each 
year 

Therefore, it is obviously important and urgent to make a 
detailed investigation of the performance of such future detec- 
tors. In this paper we focus on how far we can probe the cos- 
mic supernova rate by SRN observations at Gd-SK and at the 
hypothetical Gd-loaded Hyper-Kamiokande (Gd-HK) detec- 
tor or Gd-loaded Underground Nucleon Decay and Neutrino 



Observatory (Gd-UNO). Because the expected event rate of 
SRNs is about 1-2 yr"' in the detectable energy range (10-30 
MeV) using a detector with the size of SK, it would be quite 
difficult to obtain the spectral information of SRNs, even if 
we observed for 5 years. On the other hand, with the currently 
proposed megaton-class detectors such as HK and UNO, we 
can expect to obtain a great deal of information about the SRN 
spectrum, which will be useful for inferring the SFR-z rela- 
tion. Using the Monte Carlo (MC) method, we simulate an 
expected SRN signal at these future detectors, and then we 
analyze these hypothetical data with a few free parameters 
and discuss implications from future SRN observations. 

This paper is organized as follows. In § |2l we give the 
formulation for calculating the SRN flux and discuss several 
models that are adopted in our calculations, and in § |3l we 
show the results of our calculation with some reference mod- 
els. In §|4]the MC simulation of the expected signal at the 
future Gd-loaded detectors, which is generated from the ref- 
erence models, is presented, and then we analyze these hy- 
pothetical data using several free parameters concerning the 
cosmic SFR. Finally, we discuss other possibilities in §0 

2. FORMULATION AND MODELS 

2.1. Formulation 

The present number density of SRN (i?,,), whose energy is 
in the interval Ei, ^ Ei, + dEi,, emitted in the redshift interval 
z ~ z+ ^^z, is given by 

dnAE,) = R5M^ +zf^dz^^^0^dEl{\ 
dz dE, 



dt dN^(E') 
= Rs^iz)—dz (1 +z)dE„ 

dz dE , 



(1) 



where E^ = {I +z)Ei, is the energy of neutrinos at redshift z, 
which is now observed as E^; Rsn(z) represents the supernova 
rate per comoving volume at z, and hence the factor (1 
should be multiplied to obtain the rate per physical volume 
at that time; dN„/dEi, is the number spectrum of neutrinos 
emitted by one supernova explosion; and the factor (1 +z)~^ 
comes from the expansion of the universe. The Friedmann 
equation gives the relation between t and z as 

dz 



-=-Hq{1+z) \/(l + n„z)(l+z)2-nA(2z+z2), (2) 
dt 

and we adopt the standard ACDM cosmology (il„, = 0.3, 51a = 
0.7, and Ho = 70 hjo km s"' Mpc"').' We now obtain the dif- 
ferential number flux of SRNs, dF^/dE^, using the relation 
dF^/dEi, = cdrin/dE^: 

^ = ^ R^^iz) '^^-^^'-^ 
dEu Ho Jo ^ dEl 



dz 



v/(l + a„z)(l+z)2-OA(2z+z2) 



(3) 



where we assume that gravitational collapses began at the red- 
shift Zmax =5. 

2.2. Model for Cosmic Star Formation Rate 

' Although we use the specific cosmological model here, the SRN flux 
itself is completely independent of such cosmological parameters, as long as 
we use observationally inferred SFR m odels (see their can cellation between 
eqs. and 0), as already discussed in lAndo et alj 120031) . 
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As our reference model for the SFR, we adopt a model 
that is based on recent progressive results of rest-frame UV, 
NIR Ha, and FIR/submillimeter observations; a simple func- 
tional form for the SFR p er unit comoving volume is given as 
jPorciani & MadaiJ2nnih 



i^*(z) = 0.32/170 



exp(3.4z) 



Mq yr ' Mpc ^ 



exp(3.8z) + 45 
v/(l + a«z)(l+z)^-»A(2z+z^) 

(1 + 2)3/2 



(4) 



Figure ^ shows the S FR ^<,(z) w ith the various data points 
from rest-frame UV ("Lillv et al.' "1996*, 'Madau et al.' 199^, 
Eteideletal. 1999), Ha line ( Gallego et al. 1995; Gronwall 
[1998; Tresse & Maddox 1998), and FIR/submilHmeter 
llFlores et al. 1999: Hughes_et_al. 19_98) observations; these 
data points are not corrected for dust extinction. In the local 



-0.5 



o 

ft 



o 








1 — ^ 


— 1 



-2.5 - 



1 2 3 

Redshift z 



Fig. 1. — Cosmic star fonnation rate as a function of redshift. Data 
poin ts are given by rest-frame UV (open triangles; ILilly et alj I199{| : 
Madau et al. 1996; Steidel et al. 1999), NIR Ha (crosses; Gallego et al 1995; 
S-onwall 1998; Tresse & Maddox 1998), and FIR/submm (filled diamonds; 
Plores et aH199ftlHughes et al. 1998) observations. The solid curve repre- 



sents our reference model given by eq. j4). The standard ACDM cosmology 
is adopted (Q„, =0.3,Qa = 0.7,//o = 70 km s"' Mpc"') 



universe, all studies show that the comoving SFR monotoni- 
cally increases with z out to a redshift of at least 1 . Our refer- 
ence model (eq. 0) is consistent with mildly dust-corrected 
UV data at low redshift; on the other hand, it may underes- 
timate the results of the other wave band observations. In 
our previous paper iAndo et alJ2003 ^. we investigated the de- 
pendence on the several adopted SFR models, which were 
only different at high-redshift regions (z > 1.5); our refer- 
ence model was referred to as the "SFl" model there. We 
showed that the SRN flux at E„ > 10 MeV is highly insen- 
sitive to the difference among the SFR models (owing to the 
energy redshift, as discussed in § 13. 1> . and therefore we do 
not repeat such discussions in the present paper 

We obtain the supernova rate (.Rsn(z)) from the SFR by as- 
suming the Salpeter IMF ((j)(m) oc m"^^^) with a lower cutoff 



around 0.5Mq, and that all stars with M > SMq explode as 
core-collapse supernovae, i.e.. 



^sn(z) = 



rl25M0 



dm 4>{m) 



jj25M0 fn(j)(m) 
= 0.0122MQV*fe)- 



(5) 

The resulting local supernova rate agrees within er- 
rors with the observed value of 7?sn(0) = (1.2 ±0.4) x 
lO'^^hL yr"' Mpc"^ (e.g., |Madau, della Valle, & Panagii 
Il998al and references therein). In fact, the totally time- 
integrated neutrino spectrum from massive stars (> 30Mq) 
could be very different from the models that we use (and 
give in the next subsection), possibly because of, e.g., black 
hole formation. However, the conversion factor appearing 
in equation (|5Jl is highly insensitive to the upper limit of 
the integral in the numerator; for instance, if we change the 
upper limit in the numerator to 25M0, the factor becomes 
O.OIOMq , which is only slightly different from the value in 
equation Q- 

2.3. Neutrino Spectrum from Supernova Explosions 

For the neutrino spectrum from each supernova, we adopt 
three reference models by different groups, i.e., simulations 
by the Lawrence Livermore (LL) group (Totanietal. 199^ 
and Thompson, Burrows, & Pinto (2003, hereafter TBP), and 
the MC study of spectral formation by Keil, Raffelt, & Janka 
(2003, hereafter KR J). In this field, however, the most serious 
problem is that the recent sophisticated hydrodynamic simu- 
lations have not obtained the supernova explosion itself; the 
shock wave cannot penetrate the entire core. Therefore, many 
points still remain controversial, e.g., the average energy ratio 
among neutrinos of different flavors, or how the gravitational 
binding energy is distributed to each flavor. All these prob- 
lems are quite serious for our estimation, since the binding 
energy released as Vg changes the normalization of the SRN 
flux, and the average energy affects the SRN spectral shape. 
Thus, we believe that these three models from different groups 
will be complementary. 

The numerical simulation by the LL group JTotani et alJ 
Il998h is considered to be the most appropriate for our estima- 
tion, because it is the only model that succeeded in obtaining a 
robust explosion and in calculating the neutrino spectrum dur- 
ing the entire burst (^ 15 s). According to their calculation, 
the average energy difference between and v^, where 
represent the nonelectron-flavor neutrinos and antineutrinos, 
was rather large and the complete equipartition of the binding 
energy was realized L^^_ = Lj^^ = L^^, where L^^ represents the 
released gravitational energy as a-flavor neutrinos. The neu- 
trino spectrum obtained by their simulation is well fitted by a 
simple formula, which was originally given by KRJ as 



dEi, 



(l + /3.)'^^-L. 

r(i+/3.)£3 



E„ 



(6) 



where E^, is the average energy; the values of the fitting pa- 
rameters for the and spectrum are summarized in Table 

in 

Although the LL group succeeded in obtaining a robust ex- 
plosion, their result has recently been criticized because it 
lacked many relevant neutrino processes that are now rec- 
ognized as important. Thus, we adopt the recent result of 
another hydrodynamic simulation, the TBP one, which in- 
cluded all the relevant neutrino processes, such as neutrino 
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Table 1. Fitting Parameters for Supernova Neutrino Spectrum 



Model 


Mass 

(Mgj) 


(MeV) 


(MeV) 


/3^„ 




(ergs) (ergs) 


Reference 


LL 


20 


15.4 


21.6 


3.8 


1.8 


4.9x10^- 5.0x10^2 


1 


TBP 


11 


11.4 


14.1 


3.7 


2.2 




2 




15 


11.4 


14.1 


3.7 


2.2 




2 




20 


11.9 


14.4 


3.6 


2.2 




2 


KRJ 




15.4 


15.7 


4.2 


2.5 




3 



References. — (n lTotanietaU1998t ri^ lThompson et all2003t r3) IKeil et all200l 



bremsstrahlung and neutrino-nucleon scattering with nucleon 
recoil. Their calculation obtained no explosion, and the neu- 
trino spectrum ends at 0.25 s after core bounce. In the strict 
sense, we cannot use their result as our reference model be- 
cause the fully time-integrated neutrino spectrum is definitely 
necessary in our estimate. However, we adopt their result in 
order to confirm the effects of recent sophisticated treatments 
of neutrino processes in the supernova core on the SRN spec- 
trum. The TBP calculations include three progenitor mass 
models, i.e., 11, 15, and 2QMq\ all of these models are well 
fitted by equation (|6j, and the fitting parameters are summa- 
rized in Table [2 The average energy for both and is 
much smaller than that by the LL calculation. Although we 
do not show this in Tabled it was also found that at least for 
the early phase of the core-collapse, the complete equiparti- 
tion of the gravitational binding energy for each flavor was 
not realized. However, it is quite unknown whether these 
trends hold during the entire burst. In this study, we adopt 
the average energy given in Tabled as our reference model, 
while we assume perfect equipartition between flavors, i.e., 
Lp^=L^_ =5.0x 10^2 ergs. 

In addition, we also use the model by KRJ. Their calcu- 
lation did not couple with the hydrodynamics, but it focused 
on the spectral formation of neutrinos of each flavor using 
an MC simulation. Therefore, the static model was assumed 
as a background of neutrino radiation, and we use their "ac- 
cretion phase model II," in which the neutrino transfer was 
solved in the background of a 150 ms postbounce model by 
way of a general relativistic simulation. The fitting parameters 
for their MC simulation is also summarized in Tabled Un- 
like the previous two calculations, their result clearly shows 
that the average energy of Vy; is very close to that of Ve. It 
also indicates that the equipartition among each flavor was 
not realized, but rather L^^ ~ L^^ ~ 2Ly^. However also in 
this case, since the totally time-integrated neutrino flux is un- 
known from such temporary information, we assume perfect 
equipartition, Lp^ = L^^ = 5.0 x 10^^ ergs, as well as that the 
average energies are the same as those in Tabled 

2.4. Neutrino Spectrum after Neutrino Oscillation 

The original Ve spectrum is different from what we observe 
as Ve at Earth, owing to the effect of neutrino oscillation. 
Since the specific flavor neutrinos are not mass eigenstates, 
they mix with other flavor neutrinos during their propagation. 
The behavior of flavor conversion inside the supernova enve- 
lope is well understood, because the relevant mixing angles 
and mass square differences are fairly well determined by 
recent solar, atmospheric, and reactor neutrino experiments. 
The remaining ambiguities concerning the neutrino oscilla- 
tion parameters are the value of ^n, which is only weakly con- 



strained (sin^ 6*13 < 0. 1 :( Apollonio et alJl999l) . and the type of 
mass hierarchy, i.e., normal (mi <C m^) or inverted {m\ ^ m^). 
We first discuss the case of normal mass hierarchy as our stan- 
dard model; in this case, the value of 6'i3 is irre levan t. The 
case of inverted mass hierarchy is addressed in § 15.31 In ad- 
dition, other exotic mechanisms, such as resonant spin-flavor 
conversion (see Ando & Sato 2003b and references therein) 
and neutrino decay (Ando 2003), which possibly change the 
SRN flux and spectrum, might work in reality. However, we 
do not consider such possibilities in this study. 

The produced P^. at the supernova core are coincident with 
the lightest mass eigenstate vi owing to the large matter poten- 
tials. Since this state vi is the lightest also in vacuum, there 
are no resonance regions in which one mass eigenstate can 
change into another state, and therefore at production ar- 
rives at the stellar surface as vi . Thus, the spectrum ob- 
served by the distant detector is 



dN^ 



We 



■+\Ue2\ -TT: ^\Ue: 



= \Ue 



, dNl 
dEy, 



-+ i-ia 



dN^ 
' ' dE,, ' 



(7) 



where the quantities with superscript represent those at 
production, Uai is the mixing matrix element between the 
a-flavor state and ith mass eigenstate, and observationally 
\Uei P = 0.7. In other words, 70% of the original survives; 
on the other hand, the remaining 30% comes from the other 
component I'x- Therefore, both the original and spectra 
are necessary for the estimation of the SRN flux and spec- 
trum; since the original spectrum is generally harder than 
that of the original D^, as shown in Tabled the flavor mixing 
is expected to harden the detected SRN spectrum. 

3. FLUX AND EVENT RATE OF SUPERNOVA RELIC NEUTRINOS 

3.1. Flux of Supernova Relic Neutrinos 

The SRN flux can be calculated by equation Q with our 
reference models given in § |2] Figure |2] shows the SRN 
flux as a function of neutrino energy for the three super- 
nova models, LL, TBP, and KRJ. The flux of atmospheric 
i>e, which becomes background events for SRN detection, 
is shown in the same fi gure lOaisser, Stanev, & Ban 198^ 
iBarr. Gaisser. & Stanev 1989). The SRN flux peaks at < 5 
MeV, and around this peak, the TBP model gives the largest 
SRN flux because the average energy of the original i>e is con- 
siderably smaller than in the other two models but the total 
luminosity is assumed to be the same. On the other hand, 
the model gives a smaller contribution at high-energy regions. 
El, > 10 MeV. In contrast, the high-energy tail of the SRN flux 
with the LL model extends farther than with the other models. 
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Fig. 2. — SRN flux calculated with three reference models of original 
neutrino spectrum: LL i.Totani et aL 1998). TBP, and KR J. The flux of at- 
mospheric neutrinos <Gaisser et alll98S ^ IBarr et alll989l) is also shown for 
comparison. 



Fig. 3. — SRN flux from various redshift ranges. LL is adopted as the 
supernova model. 



and it gives flux more than 1 mag larger at E^, = 60 MeV. This 
is because the high-energy tail was mainly contributed by the 
harder component of the original neutrino spectrum; in the 
case of the LL calculation, the average energy of the harder 
component ly^ is significantly larger than that of the other two 
calculations, as shown in Tabled We show the values of the 
SRN flux integrated over the various energy ranges in Table|2] 

(In the following, we refer only to the upper pa rt of Table 
121 the values in the lower part are discussed in § 15.31 ) The 
total flux is expected to be 1 1-16 cm"^ s"' for our reference 
models, although this value is quite sensitive to the shape of 
the assumed SFR, especially at high-z. The energy range in 
which we are more interested is high-energy regions such as 
> 19.3 MeV and E^y > 11.3 MeV, because as discussed 
below, the background events are less critical and the reaction 
cross section increases as oc /s^- In such a range, the SRN flux 
is found to be 1.3-2.3 cm'^ s"' (E^ > 11.3 MeV) and 0.14- 
0.46 cm"^ s"' (El, > 19.3 MeV). Thus, the uncertainty about 
the supernova neutrino spectrum and its luminosity gives at 
least a factor 2-4 ambiguity to the expected SRN flux in the 
energy region of our interest. 

Figure shows the contribution by supernova neutrinos 
emitted from various redshift ranges. At high-energy region 
El, > 10 MeV, the dominant flux comes from the local super- 
novae (0 < z < 1), while the low-energy side is mainly con- 
tributed by the high-redshift events (z > 1). This is because 
the energy of neutrinos that were emitted from a supernova 
at redshift z is reduced by a factor of (1 +z)"' reflecting the 
expansion of the universe, and therefore high-redshift super- 
novae only contribute to low-energy flux. We also show the 
energy-integrated flux from each redshift range in Table |2l in 
the case of the LL supernova model. From the table, it is 
found that in the energy range of our interest, more than 70% 
of the flux comes from local supernova explosions at z < 1 , 



while the high-redshift (z > 2) supernova contribution is very 
small. 

3.2. Event Rate at Water Cerenkov Detectors 

The water Cerenkov neutrino detectors have greatly suc- 
ceeded in probing the properties of neutrinos as elementary 
particles, such as neutrino oscillation. The SK detector is one 
of these detectors, and its large fiducial volume (22.5 kton) 
might enable us to detect the diffuse background of SRNs. 
Furthermore, much larger water Cerenkov detectors such as 
HK and UNO are being planned. SRN detection is most 
likely with the inverse /?-decay reaction with protons in wa- 
ter, P ffP — > e'^n, a nd its cross se ction is precisely understood 
( Vogel & Beacom 1999; Strum ia & Vissa ni 2003). In our cal- 
culation, we use the trigger threshold of SK-I (before the ac- 
cident). 

The expected event rates at such detectors are shown in Fig- 
ures |4| and |5] in units of (22.5 kton yr)"' MeV"'; with SK, it 
takes a year to obtain the shown SRN spectrum, while with 
HK and UNO, much less time [1 yr x (22.5 kton/Vfid), where 
Vfid is the fiducial volume of HK or UNO] is necessary be- 
cause of their larger fiducial volume. Figure compares 
the three models of the original supernova neutrino spectrum, 
and Figure [S] shows the contribution to the total event rate 
from each redshift range. In Table|2]we summarize the event 
rate integrated over various energy ranges for three supernova 
models. The expected event rate is 0.97-2.3 (22.5 kton yr)"' 
for Ee > 10 MeV and 0.25-1.0 (22.5 kton yr)"' for E, > 18 
MeV. This clearly indicates that if the background events that 
hinder the detection are negligible, the SK has already reached 
the required sensitivity for detecting SRNs; with the future 
HK and UNO, a statistically significant discussion would be 
possible. This also shows that the current shortage of our 
knowledge concerning the original supernova neutrino spec- 
trum and luminosity gives at least a factor of 2 (E^, > 10 MeV) 
to 4 {E,y > 1 8 MeV) uncertainty to the event rate at the high- 
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Table 2. Flux and Event Rate of Supernova Relic Neutrinos 



Flux [cm ^ s '] Event Rate [(22.5 kton yr) '] 



Model Redshift Range Total £^>11.3MeV £^ > 19.3 MeV > 10 MeV £(,>18MeV 



LL Total 11.7 2.3 0.46 2.3 1.0 

0<z<l'' 4.1 (35.3) 1.6(70.9) 0.39 (85.2) 1.7(77.5) 0.9 (87.5) 

l<z<2= 4.9(42.0) 0.6(26.3) 0.06(14.0) 0.5(20.6) 0.1 (11.9) 

2<z<3'' 1.8(15.1) 0.1(2.5) 0.0(0.7) 0.0(1.7) 0.0(0.5) 

3<z<4'' 0.6(5.3) 0.0(0.2) 0.0(0.0) 0.0(0.1) 0.0(0.0) 

4<z<5'' 0.2(2.1) 0.0(0.0) 0.0(0.0) 0.0(0.0) 0.0(0.0) 

TBP Total 16.1 1.3 0.14 0.97 0.25 

KRJ Total 12.7 2.0 0.28 1.7 0.53 



Inverted Mass Hierarchy with Large 

Tl Total 9A Jl 094 11 23 

TBP Total 13.8 1.9 0.30 1.6 0.58 

KRJ Total 12.4 2.2 0.38 2.0 0.76 



Note. — Values in the upper part are evaluated for the case of normal mass hierarchy (or inverted mass hierarchy with sufficiently small 6u, i.e., sin 2613 < 10"'), which we use as 
our standard model. On the other hand, values in the lower part are applicable only when the value of 6^13 is large enough to induce completely adiabatic resonance, i.e., sin" 2^13 ^ 10^, 
in the case of inverted mass hierarchy. 

"Contributions from each redshift range to the total (0 < z < 5) value are shown in parentheses as percentages. 




Positron Energy [MeV] 



Fig. 4. — Event rate at water Cerenkov detectors in units of (22.5 kton yr) 
for three supernova models. 



energy range (actual detection range). We also summarize the 
contribution from each redshift range in the same table, es- 
pecially for the calculation with the LL model. The bulk of 
the detected events will come from the local universe (z < 1), 
but the considerable flux is potentially attributed to the range 
l<z<2. 

3.3. Comparison with Other Studies and Current 
Observational Limits 

There are many past theoretical researches concerning the 
SRN flux estimation based on a theoretically/observati onally 
modeled cosmic SFR jTotani et alJ Il996t iMalanevI 119971 




10 20 30 40 50 60 



Positron Energy [MeV] 

Fig . 5 . — Event rate at water Cerenkov detectors in units of (22.5 kton yr) 
from various redshift ranges. LL is adopted as the supernova model. 



iHartmann & WooslevI fT997t lAndo et alJ l200l . Here we 
briefly compare our results obtained in § 13. H and § 13. 21 with 
these past analyses. Our basic approach in the present pa- 
per is the same as that in Ando et al. (2003), in which the 
LL supernova model was adopted. Thus the values for the 
L L model given in Table |2l are almost the same as those 
inlAndo et all j2003h Two other studies JTotani etalJI 19961: 
IHartmann & Wooslevlll99'7h also used a similar SFR-z rela- 
tion at low-redshift, and therefore their results are very con- 
sistent with the present one (the LL model) at high-energy re- 
gion > 10 MeV. Since the SFR model adopted bv MalanevI 
i 1997b gave a rather lower value at low-redshift, the result- 



STAR FORMATION RATE AND SUPERNOVA NEUTRINOS 



7 



ing SRN flux at high-energy regions was about a factor 2 
smaller than our LL model or the other on es ( Totani et al. 
Il996t iHartmann & WoosIevlfT997t lAndoetaL.2003.) . Thus, 
our calculation with the LL supernova model gives values 
quite consistent with past studies within a factor of 2, but all 
of those studies used the supernova model that is very sim- 
ilar to the LL model. In fact, the present study is the first 
one to investigate the dependence on the adopted supernova 
models, by using the various original spectra of the different 
groups (LL, TBP, and KRJ). As already mentioned in § 13.11 
it was found that the ambiguity concerning the original neu- 
trino spectrum varies the resulting value of the flux by at least 
a factor of 2. 

In additi on, there are severa l other studies o n the 
SRN flux (Totani & Satd [T99l 'Kanlinshat et al.' lood; 
Eukugita & Kawasaki 200%. As for Totani & Sato ( 1995), 
the authors used a constant supernova rate model in order 
to investigate the dependence on cosmological parameters; 
they gave a very large value (~ 3 cm"^ s~' at Ei, > 19.3 
MeV), which is already excluded observationally, because 
they adopted a rather large supernova rate. Concerning the 
other two studies, since neither of them gives a specific value 
for the SRN flux, we cannot compare ours with theirs; they 
focus ed on giving theoretical upper limit ( Kaplinahat et al. 
I2000h or probing the cosmic SFR with the current observa- 
tional upper limit by SK (Fukugita & Kawasaki 2003). 

Observationally, the SK collaboration gave a very stringent 
upper hmit to the SRN flux at > 19.3 MeV, i.e., 1.2 cm"^ 
s"' (90% CL; Malek et al. 200S- This number can be directly 
compared with our predictions summarized in Table |2l Our 
predicted values are 0.46, 0.14, and 0.28 cm"- s"' for the LL, 
TBP, and KRJ models, respectively. Thus, the current SK up- 
per limit is about a factor 2.5-8.5 larger than our predictions 
with the reference model for the cosmic SFR, depending on 
the adopted original neutrino spectrum. 

3.4. Background Events against Detection 

In § 13.21 we calculated the expected SRN spectrum at the 
water Cerenkov detectors on Earth, but the actual detection is 
quite restricted because of the presence of other background 
events. There are atmospheric and solar neutrinos, antineu- 
trinos from nuclear reactors, spallation products induced by 
cosmic -ray muons, and decay products of invisible mu ons (for 
a deta iled discussion of these backgrounds, see An do et al.l 
I2003h . For the pure-water Cerenkov detectors, there is no en- 
ergy window in which the flux of any backgrounds is much 
smaller than the SRN flux. 

However, as proposed by iBeacom & VaginsI J2003I) . if we 
use Gd-loaded detectors, the range 10-30 MeV would be an 
energy window because we can positively distinguish the 
signal from other backgrounds such as solar neutrinos (i^^), in- 
visible muon events, and spallation products; this is realized 
by capturing neutrons that are produced by the D^p interac- 
tions. Above 30 MeV, the SRN flux becomes smaller than the 
flux of atmospheric P^, as shown in Figure |2l because they 
are of the same flavor, it is in principle impossible to distin- 
guish them from the SRN D,,. On the other hand, below 10 
MeV the reactor neutrinos (Pg) are the dominant background 
in the case of SK or HK; because the flux of reactor neutrinos 
strongly depends on the detector site, it may be possible to 
further reduce this lower energy cutoff (10 MeV) in the case 
of UNO. 

The neutron capture efficiency by Gd is estimated to be 
90% with the proposed 0.2% mixture by mass of GdCh in 



water (iBeacom & VaginsI 120031) . and subsequently 8 MeV 
gamma cascade occurs from the excited Gd. The single- 
electron energy equivalent to this cascade was found to be 3- 
8 MeV by careful simulation (Hargrove et al. 1995), and with 
the trigger threshold adopted in SK-I, only about 50% of such 
cascades can be detected actually. However, it is expected that 
SK-III, which will begin operation in mid-2006, will trigger 
at 100% efficiency above 3 MeV, with good trigger efficiency 
down to 2.5 MeV (Beacom & Vagins 2003). In that case most 
of the gamma cascades from Gd will be detected with their 
preceding signal of positrons. From this point on, we assume 
100% efficiency; even if we abandon this assumption, it does 
not affect our physical conclusion, since the relevant quantity 
representing the detector performance is (fiducial volume) x 
(time) X (efficiency), which we call effective volume. 

4. MONTE CARLO SIMULATION FOR FUTURE DETECTOR 
PERFORMANCE 

In this section we predict the expected signal at future de- 
tectors, such as Gd-SK, Gd-HK, and Gd-UNO, using the MC 
method with our reference models. These pseudodata are then 
analyzed using several free parameters concerning the SFR. 
Although we focus on how far the SFR can be probed by 
SRN observation, the uncertainty from the supernova neu- 
trino spectrum would give a fair amount of error. However, 
this problem can be solved if a supernova explosion occurs in 
our Galaxy; the expected event number is about 5000-10,000 
at SK, when supernova neutrino burst occurs at 10 kpc, and it 
will enable a statistically significant discussion concerning the 
neutrino spectrum from supernova explosions. Even if there 
are no Galactic supernovae in the near future, remarkable de- 
velopment of the supernova simulation can be expected with 
the growth of computational resources and numerical tech- 
nique. With such developments, the supernova neutrino spec- 
trum and luminosity may be uncovered, and the ambiguity is 
expected to be reduced significantly. Thus, in this paper we 
assume that the supernova neutrino spectrum is well under- 
stood and that our reference models are fairly good represen- 
tatives of nature; we analyze the SFR alone with several free 
parameters. 

The basic procedure of our method is as follows. (1) We 
simulate the expected signal (spectrum) at a Gd-loaded detec- 
tor in the range 10-30 MeV, assuming that there are no back- 
ground events. In that process, we use our reference models 
for the generation of the SRN signal (eq. for the SFR and 
the LL model as neutrino spectrum). (2) Then we analyze the 
SRN spectrum using the maximum likelihood method with 
two free parameters of the SFR and obtain a set of the best- 
fit values for those parameters; they are concerned with the 
supernova rate as 



^sn(z) ■■ 



ict dO 



for z < 1 , 
for z > 1 , 



(8) 



where represents the local supernova rate and a deter- 
mines the slope of supernova rate evolution. Although it is 
recognized that the SFR-z relation increases from z = to z = 1 
from various observations using light, the actual numbers for 
the absolute value and the slope of the SFR-z relation are still 
a matter of controversy and independent confirmation, such 
as ours, is needed. We assume that the comoving SFR is con- 
stant at z > 1; even if we changed this assumption, the re- 
sult would be the same because the bulk of the detected event 
comes from local supernovae. (3) We perform 10^ such MC 
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simulations and obtain 10^ independent sets of best-fit param- 
eters. Then we discuss the standard deviation of the distribu- 
tions of such best-fit parameter sets and the implications for 
the cosmic SFR. 

4.1. Performance of the Gd-loaded Super-Kamiokande 
Detector 

In this subsection we discuss the performance of Gd-SK 
for 5 years, or an effective volume of 22.5 kton x 5 yr. Be- 
cause the expected event number is only 10, the parameters 
rIj^ and a cannot both be well determined at once. There- 
fore, we fix one of those parameters with some value inferred 
from other observations. First, the value of R^j^ was fixed to 
be 1.2 X 10""* yr"' Mpc~', which was inferred from the lo- 
cal supernova survey (Madauet al. 1998a), and we obtained 
the distribution of the best-fit values of parameter a. Figure 
|6l shows the expected SRN spectrum; points with error bars 
represent the result of one MC simulation, and the dashed his- 
togram is the spectrum with the best-fit parameter {a = 2.97). 
Thus, from one realization of the MC simulation, we obtain 
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Fig. 6. — Expected SRN spectrum at a detector of effective volume 22.5 
kton X 5 yr The data points represent tlie result of MC simulation, and the 
error bars include statistical errors alone. These data generated by MC sim- 
ulation were analyzed assuming = 1.2 X IQ-'' yr-' Mpc"', and using a 
as a free parameter The best-fit value for a is 2.97 and it resulted in dashed 
histogram in this figure. 



one best-fit parameter The result of lO"' MC simulations are 
shown in Figure0as a histogram of the distribution of best-fit 
parameters a {solid histogram). The average value of these 
10"* values for a is found to be 2.67, and the standard devi- 
ation is 0.80, i.e., a = 2.67 ±0.80. A no-evolution (constant 
supernova rate) model would be excluded at the 3.3(7 level 
from the SRN observation alone with an effective volume of 
22.5 kton x 5 yr. 

Then in turn, we fixed the value of a to be 2.9 in order to 
obtain the distribution of best-fit values for the local super- 
nova rate from the SRN observation. The result of 10^ 
MC generations and analyses in this case is shown in Fig- 
ure |8] The average value for /Jg^^ is 1.2 x lO"'' yr"' Mpc"^, 



Fig. 7. — Distribution of lO' best-fit values for a, which are obtained from 
the analyses of each MC generation. The effective volume is 22.5 kton X 5 
yr for solid histogram, and 440 kton X 5 yr for dashed histogram. The value 



of the local supernova rate is fixed to be Ri^ = 1.2 X 10 yr ' Mpc ^ 
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Fig. 8. — Distribution of 10^ best-fit values for local supernova rate Rgp^, 
which are obtained from the analyses of each MC generation. The effec- 
tive volume is 22.5 kton X 5 yr for solid histogram, and 440 kton X 5 yr for 
dashed histogram. The value of a is fixed to be 2.9. 



^SN 



and the standard deviation is 0.4 x 10 yr Mpc , i.e., 
(1.2±0.4)x 10-"* yr"' Mpc'l 
The results obtained with the above calculations are sum- 
marized in Table 13 In Figure |5] we compare the super- 
nova rate model in which the parameter is inferred from the 
MC simulations with the "true" reference model; the cases of 
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Table 3. 
Model 



Expected Sensitivity of Future Detectors to Supernova Rate 



Detector 


Effective Volume 
(22.5 kton yr) 


Fixed 
Parameter 


a 


5a /{a) 
(%) 


(I0-* yr-' Mpc-3) 


(%) 


SK 


5 




2.7 ±0.8 


30.0 


1.2 (fixed) 






5 


a 


2.9 (fixed) 




1.2 ±0.4 


28.3 


HK or UNO 


97.8 




2.5 ±0.2 


7.8 


1.2 (fixed) 






97.8 


a 


2.9 (fixed) 




1.0±0.1 


7.7 




97.8 




3.5± 1.3 


36.7 


0.88 ±0.48 


54.8 



fixed TJgj^ and a are shown in Figures|9j? andl^J?, respectively. 
The allowed region at the la level is located between the two 



o 



10- 



10- 



03 
K 

> 
O 

a 

CD 

a 

CO 



10- 



10- 



1 ' 1 ' 
(a) Fixed 


1 






1 1 1 1 1 , 

(b) Fixed a 


^ \ ^ 


- , 1 , 1 , 


1 , - 



0.2 



0.4 



0.6 O.i 



Redshift z 



Fig. 9. — Supernova rate as a function of redshift. In both panels, solid 
curves represent our reference model, (a) The allowed region at the la level, 
concerning the fitting parameter a with fixed Rg^^, is shown as the area be- 
tween the two dotted curves for an effective volume of 22.5 kton X 5 yr 
and as the area between the two dashed curves for an effective volume of 
440 kton X 5 yr. (b) Same as (a) but for fitting parameter R'^^ with fixed a. 



dotted curves, while the solid curve represents our reference 
model. Thus, with the Gd-SK detector we can roughly repro- 
duce the supernova rate profile at z < 1 for 5 years operation, 
although it is still statistically insufficient. 

4.2. Future Gd-loaded Megaton-Class Detectors 

We consider future megaton-class detectors such as Gd-HK 
or Gd-UNO. With these detectors, the effective volume that 
we consider, 440 kton x 5 yr, is expected to be realized in 
several years from the start of their operation. First we did 
the same analysis adopted in the previous subsection, i.e., we 
fixed one of relevant parameters, a or R\^, and investigated 
the dependence on the remaining parameter The values that 
we used for fixed parameters were the same as those given 
in the previous subsection. The result of these cases are also 
shown in Figures [7] and |S] as dashed histograms, which give 
a = 2.51±0.20and7;gN = (l-0±0.1)x lO""* yr"' Mpc-\ re- 
spectively, and these values are also summarized in Table |3l 



The statistical errors are considerably reduced compared with 
the case of 22.5 kton x 5 yr, because of the ^ 20 times larger 
effective volume. Thus, future magaton detectors will possi- 
bly pin down, within 10% statistical error, either the index of 
supernova rate evolution a or the local supernova rate R^^ if 
the other is known in advance. The dashed curves in Figure|9] 
set the allowed region of the supernova rate at the 1 a level 
by the considered detectors, well reproducing the assumed 
model. 

In principle, we can determine both parameters by SRN ob- 
servation, because R^^ is concerned with the absolute value 
of the flux alone but a is concerned with both the absolute 
value and the spectral shape; i.e., these two parameters are 
not degenerate with each other. Thus, we repeated the same 
procedure but without fixing the values of a or R^^. The dis- 
tribution of lO-' best-fit parameter sets of is shown 
in Figure [TO] for a detector with an effective volume of 440 
kton X 5 yr; the mean values and the standard deviations are 
a = 3.5 ± 1.3 and 7?°^ = (8.8 ±4.8) x 10"^ yr"' Mpc'l Even 




Fig. 10. — Distribution of 10 best-fit parameter sets (o,/?gj^). Each dot 
represents the result of one MC generation and an accompanying analysis. 
The effective volume is 440 kton X 5 yr 



though the effective volume is as large as 440 kton x 5 yr, it is 
still insufficient for determining both parameters at once. For 
another trial, we also carried out the same MC simulations, 
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but using a hypothetical (and unrealistic) effective volume as 
large as 440 kton x 10"* yr In that case the free parameters 
are found to be quite well constrained at a = 3.68 ± 0.03 and 
RIj^ = (7.1 1 ± 0.09) X 10"^ yr-i Mpc'^ 

5. DISCUSSION 

5.1. Supernova Rate at High-Redshift Region 

At the detection energy range 10-30 MeV that we have con- 
sidered, the main contribution to the SRN event rate comes 
from low-redshift region < z < 1, as shown in Figure|5]and 
Table |2] However, if we can reduce the lower energy thresh- 
old Eth, we expect that the contribution of supernova neutrinos 
from high-redshift z > 1 becomes enhanced. The value of Eth 
is restricted to 10 MeV because at energy regions lower than 
this, there is a large background of reactor neutrinos; its re- 
moval is impossible with the current detection methods. Since 
the SK and HK detectors are and will be located at Kamioka 
in Japan, they are seriously affected by background neutrinos 
from many nuclear reactors. If some large-volume detectors 
were built at a location free from such background, the lower 
threshold energy could be reduced, enabling us to probe the 
high-redshift supernova rate. In this subsection, thus, we dis- 
cuss the detector performance as a function of the value of 
Eth- 

In Figure [TTl; we show three toy models of comoving den- 
sity of supernova rate as a function of redshift. These models 
exactly coincide with each other at z < 1 (and also with the 
previous reference model represented by eq. 0) but seri- 
ously differ at z > 1 . We calculate the expected event number 




Threshold Energy E^^ [MeV] 



Fig. 1 1 . — (a) Three toy models for comoving density of supernova rate as 
a function of redshift. (b) Expected event number A' at Eth < Ee < 30 MeV 
as a function of £th. for the fiducial volume of 440 kton X 5 yr The line types 
correspend to those in (a). The upper and lower curves of each type represent 
N+\/N and N—\/N, respectively; i.e., the area between the two curves is the 
allowed region at the 1 a level. 

for Eth <Ec <3Q MeV at a detector with an effective volume 
of 440 kton x 5 yr, using these toy models, the LL spectrum, 
and trigger threshold expected at SK-111. The result is shown 
in FigureFTTt'. As expected, the discrepancy among the three 



models becomes larger as we reduce the threshold energy. In 
particular, the model that produces larger numbers of super- 
novae at z > 1 (solid curve) is satisfactorily distinguishable in 
the case of sufficiently low E±. This is because the larger su- 
pernova rate at high-redshift region z > 1 increases the frac- 
tion of its contribution to the SRN flux. On the other hand, 
the model with a lower supernova rate relatively increases 
the contribution from low-redshift region z < 1, and therefore 
the difference between constant (dotted curve) and decreasing 
model (dashed curve) is less prominent. 

5.2. Probing Supernova Neutrino Properties 

Until this point, we have assumed that the properties of su- 
pernova neutrinos, such as the average energy difference be- 
tween flavors and luminosities of neutrinos of different fla- 
vors, will be quite well understood when future SRN detection 
comes within reach. However, this assumption itself is quite 
unclear because Galactic supernova explosions, which would 
give us rich information on the supernova neutrino spectrum 
and luminosity, may not occur by the time we are ready for 
the SRN detection. Furthermore, there is no assurance that 
the numerical experiments will succeed in obtaining the su- 
pernova explosion itself and predicting the supernova neutrino 
properties precisely by then. Thus, SRN observation might be 
the only probe of the supernova neutrino properties. 

In this subsection, we discuss how far we can derive the su- 
pernova neutrino properties from SRN observation. We have 
already shown that even using data of 440 kton x 5 yr, at most 
only two free parameters can be satisfactorily constrained. 
Therefore, we now have to adopt another assumption such that 
the evolution of the supernova rate is quite well understood 
by future observations with the various planned satellites and 
telescopes. The procedure is basically the same as that of the 
previous section; i.e., we run 10^ MC simulations and analyze 
these psudodata to obtain the best-fit values for two free pa- 
rameters, £p and EyjEjj^,. The values of [3y and defined 
in equation |6j are assumed to be = 4.0, (3^, = 2.2, and 
Lp^ =Lix^ = 5.0 x 10''~ergs. As a result of such calculations, we 
obtain the distribution of the two parameters, which is charac- 
terized by = (15.9 ± 1 .3) MeV and E^JE;,^ = 1 .5 ±0.4; al- 
though this well reproduces the LL model, the errors are still 
very large. Considering that many uncertainties concerning 
the SFR estimate possibly remain even in future updated ob- 
servations, the errors to these quantities would be much larger 
than the purely statistical ones given above. 

5.3. Inverted Mass Hierarchy 

Throughout the above discussion, we have assumed normal 
hierarchy of neutrino masses (mi <C mj). However, the case 
of inverted mass hierarchy has not been experimentally ex- 
cluded yet, and we explore this possibility in this subsection. 
In this case, flavor conversions inside the supernova envelope 
change dramatically, com pared with the normal mass hierar- 
chy already discussed in § 12.41 Since is the lightest, are 
created as P3, owing to large matter potential. In that case, it 
is well known that at a so-called resonance point, there occurs 
a level crossing between i/| and Dj, (for a more detailed dis- 
cussion, see, e.g..lDi ghe & Smir nov 2000). At this resonance 
point, complete Di ^ P3 conversion occurs when the so-called 
adiabaticity parameter is sufficiently small compared to unity 
(it is said that resonance is "nonadiabatic"), while conversion 
never occurs when it is large (adiabatic resonance). The adia- 
baticity parameter 7 is quite sensitive to the value of 613, i.e.. 
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7 oc sin^26'i3; when sin^26li3 > 10"^ (sin226'i3 < 10"^), the 
resonance is known to be completely adiabatic (nonadiabatic) 
dPighe & Smirnovl2000h . When the resonance is completely 
nonadiabatic (because of small ^o), the situation is the same 
as in the case of normal mass hierarchy akeady discussed in 
§ 12.41 (because at production become v\ at the stellar sur- 
face), and the spectrum after oscillation is represented by 
equation Q. On the other hand, adiabatic resonance (due to 
large forces at production to become P3 when they es- 
cape from the stellar surface, and therefore the observed 
spectrum is given by 



d^Ej _ dill ^ 



(9) 



The second equaUty follows from the fact that the value of 
|f/e3p is constrained to be much smaller than unity from re- 
actor experiments ( Apollonio et al. 1999). Thus, equation (|9} 
indicates that complete conversion takes place between Vg and 
Vx- When the value of ^13 is large enough to induce adiabatic 
resonance (sin^20i3 > 10~^), the obtained SRN flux a nd sp ec- 
trum should be very different from ones obtained in §§ 13. H and 
13.21 The SRN flux and event rate in this case were calculated 
with equations Q and and the results are summarized in 
the lower part of Table |2 The values (with the LL model) 
shown i n this table are consis tent with the previous calcula- 
tion by Ando & Sato (2003a), in which numerically calcu- 
lated conversion probabilities were adopted with some spe- 
cific oscillation parameter sets (which include a model with 
inverted mass hierarchy and sin^20i3 = 0.04), as well as real- 
istic stellar density profiles. 

The total flux becomes 9.4-14 cm"^ s"', somewhat smaller 
than the values given in the upper part of the same table, be- 
cause the total flux is dominated by the low-energy region. 
The fluxes at > 19.3 MeV are enhanced to be 0.30-0.94 
cm"^ s"' , but this is still below the current 90% CL upper limit 
of 1 .2 cm"^ s"' obtained by the SK observation. The event 
rate at the future detectable energy range, E^ > 10 MeV, is 
expected to become 1.6-3.8 yr~', which is considerably larger 
than the values in the case of normal mass hierarchy, 0.97-2.3 
yr"' . The increase (decrease) of the flux and event rate inte- 
grated over the high (total) energy region is due to the very 
high efficiency of the flavor conversion, — *■ D^, inside the 
supernova envelope; because the original are expected to 
be produced with larger average energy as shown in Table 
d] the efficient conversion makes the SRN spectrum harder, 
which enhances the flux and event rate at the high-energy re- 
gion. Thus, if the inverted mass hierarchy, as well as the large 
value for 6*13, were realized in nature, SRN detection would 
be rather easier, compared with the other cases. Although we 
do not repeat the MC simulations that were introduced in §|3 
the results can be easily inferred; the statistical errors in this 
case would be ^ (3.8/2.3)'''^ =1.3 times smaller than the val- 
ues given in Table |3l because they are inversely proportional 
to the square root of the event number. 

6. CONCLUSIONS 

In the present paper, we have investigated the flux and event 
rate of SRNs and discussed their implications for the cos- 
mic SFR. Since SRNs are diffuse neutrino background emit- 
ted from past core-collapse supernova explosions, they con- 
tain fruitful information not only on the supernova neutrino 
spectrum itself but also on the supernova rate in the past and 
present universe, which is quite difficult to estimate because. 



e.g., the problem of dust extinction is nontrivial. As refer- 
ence models, we adopted the supernova rate model based on 
recent SFR observations (eq. |4J) and the supernova neutrino 
spectrum numerically calculated by three groups (LL, TBP, 
and KRJ). As a result of our calculations, the flux integrated 
over the entire energy region was found to be 12-16 cm"^ 
s"', depending on the adopted supernova neutrino spectrum 
(Table |2}. Although there is no energy window for the SRN 
detection at present owing to various background events, in 
the near future, it is expected that the energy region of 10-30 
MeV will be utilized for SRN detection. This is due to the 
technique of neutron capture by dissolved Gd. In the detec- 
tion energy range /s^ > 10 MeV, the SRN event rate was found 
to be 0.97-2.3 yr"' at a detector with a fiducial volume of 22.5 
kton (Table |2). 

We also simulated the expected signal with one set of the 
reference models by using the Monte Carlo method and then 
analyzed these pseudodata with several free parameters, ob- 
taining one set of best-fit values for them. MC simulations 
repeated 10^ times gave lO"* independent best-fit parameter 
sets, and we gave a statistical discussion using their distribu- 
tion. First of all, we used parameterization such that /?sn(z) = 
^sn(1+z)"' where and a are free parameters, assuming 
that the supernova neutrino spectrum and luminosity are well 
understood by way of a future Galactic supernova neutrino 
burst or future development of the numerical supernova sim- 
ulations. The obtained distribution for these two parameters 
was found to be represented by a = 2.7 ±0.8, Sa/ {a) = 30% 
and Rl^ = (1.2 ±0.4) x lO'^ yr"' Mpc'^SRl^/ (RI^) = 28% 
for a detector with an effective volume of 22.5 kton x 5 yr, 
and a = 2.5 ±0.2, Sa/ {a) = 7.8% and 7?'^^ = (1-0 ±0.1) x 
10"'* yr"' Mpc-\SRl^/(Rlj^) = 7.7% for a detector with an 
effective volume of 440 kton x 5 yr, where one of the param- 
eters is fixed (Figs. 0and|8j Table |3ji. The parameterized 
supernova rate models with the obtained parameter values are 
compared with the assumed reference model in Figure|9j and 
we found that the fitting model well reproduced the reference 
model. On the other hand, if we fix neither value for these 
two parameters, the expected errors become rather large at 
Sa/(a) = 37% and (5/?'^n/(^sn> = 55%, even with an effective 
volume of 440 kton x 5 yr 

In addition, we explored several other possibilities in § |5l 
First, we discussed the dependence of the event number on 
the adopted lower cutoff energy. Although below 10 MeV 
there is a background of reactor neutrinos, their flux strongly 
depends on the detector sites, and the lower energy threshold 
Eih could possibly be reduced. We investigated the expected 
event number for E^h < Eg < 30 MeV as a function of E^h in 
Figure^Jfor various toy models of supernova rate and found 
that the model that produces larger number of supernovae at 
z > 1 is satisfactorily distinguishable in the case of sufficiently 
small Eth. Second, the SRN spectrum as a potential probe of 
the supernova neutrino spectrum itself was investigated, be- 
cause such an approach might be very important if there are 
no Galactic supernova explosions in the near future or no suc- 
cessful numerical supernova simulations. We discussed us- 
ing the same MC procedure, but assuming that the supernova 
rate is quite well understood. Although the obtained distri- 
bution reproduces properties of the LL spectrum, the errors 
were still found to be large, and considering the uncertainties 
concerning the SFR, these errors are only lower limits; the 
actual errors would be much larger. Finally, the case of an in- 
verted mass hierarchy was investigated. We showed that only 
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in the case in which sin^20i3 > 10 ^ the values of the SRN 
flux should be modified. The results in the case of completely 
adiabatic resonance, which is realized when sin^lOn > lO""*, 
are shown in the lower part of Table |2] In this case, it was 
found that the expected event rate would be enhanced to 1.6- 
3.8 yr~', although these values are still below the current up- 



per bound; SRN detection would be more probable in this 
case. 
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